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Abstract 

We make a new estimate of the extragalactic radio background down to 
kHz frequencies based on the observed luminosity functions and radio spectra 
of normal galaxies and radio galaxies. We have constructed models for the 
spectra of these two classes of objects down to low frequencies based on 
observations of our Galaxy, other normal galaxies and radio galaxies. We 
check that the models and evolution of the luminosity functions give source 
counts consistent with data and calculate the radio background expected 
from kHz to GHz frequencies. 

The motivation for this calculation is that the propagation of ultra-high 
energy gamma-rays in the universe is limited by photon-photon pair produc- 
tion on the radio background. Electromagnetic cascades involving photon- 
photon pair production and subsequent synchrotron radiation in the inter- 
galactic magnetic field may develop. Such gamma-rays may be produced 
in acceleration sites of ultra-high energy cosmic rays, as a result of inter- 
actions with the microwave background, or emitted as a result of decay or 
annihilation of topological defects. We find that photon-photon pair produc- 
tion on the radio background remains the dominant attenuation process for 
gamma-rays from 3 x 10 10 GeV up to GUT scale energies. 
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1 Introduction 



The universe is not transparent to high energy gamma-rays due to inter- 
actions with low energy photons of the extragalactic radiation fields, the 
most important process being photon-photon pair production. For example, 
interactions in the cosmic microwave background radiation give a mean in- 
teraction length of less than 10 kpc at 10 6 GeV as has been known since 
soon after the discovery of the microwave background ]l|, [|. The threshold 
for interactions on the microwave background is ~ 10 5 GeV, and at lower 
energies interactions on the infrared and optical backgrounds limit the trans- 
parency at TeV energies (e.g. || [|]) Other components of the extragalactic 
background radiation are discussed in the review of Ressel and Turner ||. 

Above ~ 10 10 GeV interactions with the radio background become more 
important than the microwave background in limiting the transparency of 
the universe to gamma-rays and controlling any resulting electromagnetic 
cascades. Both the infrared and radio backgrounds are poorly known due to 
our location within our Galaxy which emits and absorbs at these wavelengths. 
The radio background was measured over twenty-five years ago 0, but the 
fraction of this radio background which is truly extragalactic, and not con- 
tamination from our own Galaxy, is still debatable. A theoretical estimate M 
was made about the same time which gave a quite different spectrum, par- 
ticularly at low frequencies. In recent cascade calculations [|], ||, [K| [ll|] the 
estimate of ref. || has been used. It is this very uncertain radio background 
which will provide target photons for UHE 7-rays above ~ 10 10 GeV. 

While gamma-ray astronomy is not currently undertaken at ~ 10 10 GeV 
energies, it is important to know the photon-photon mean interaction length 
at these energies because cascading involving gamma-rays at these energies 
occurs in top-down models for the origin of the highest energy cosmic rays. 
The highest energy cosmic rays have energies of 200 EeV O, IT3 and 300 



EeV [R], and are well above the "Greisen-Zatsepin-Kuzmin cut-off" [15, 16 
at 50 EeV in the spectrum of cosmic ray protons due to pion photoproduction 
in the microwave background, which is expected if the cosmic rays originate 
further than a few tens of Mpc (e.g., ]17|, |18[). Of the various models proposed 



to account for the origin of these high energy cosmic rays |L9[ [L8|, £0[ [21], ^] , 
one of the more tantalizing speculations is that the highest energy cosmic 
rays may be ultimately due to the decay of supermassive X particles [E3 



2J, g^, |26| , themselves radiated during collapse or annihilation of topological 
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defects, remnants of an early stage in the evolution of the Universe. The 
X particles have GUT-scale masses of ~ 10 16 GeV or lower, depending on 
the theory, and decay into leptons and quarks at lower energies. The quarks 
themselves fragment into a jet of hadrons which, it is supposed, could produce 
the highest energy cosmic rays, although there is some debate as to whether 
a sufficiently large fraction of the energy of the defect could end up in high 
energy particles [p7fl . In any case, much of the radiation is likely to emerge 
in the electromagnetic channel and initiate an electromagnetic cascade in 
the ambient radiation field, in which collisions with radio photons play an 
important role. 

In this paper we make a new calculation of the extragalactic radio back- 
ground down to kHz frequencies based on the infrared luminosity function of 
normal galaxies recently determined from IRAS source counts, the observed 
radio-infrared correlation, and the luminosity function of radio galaxies, to- 
gether with recent models for radio spectra of these objects. Finally, we 
calculate the mean free path for 7-rays in the extragalactic radio background 
radiation. 



2 Calculation of radio background 

The main contributions to the radio background will be from normal galaxies 
and radio galaxies and we will discuss the radio spectra of these objects down 
to kHz frequencies and construct models for their spectra. Using appropri- 
ate luminosity functions, we will then integrate over luminosity and redshift 
to obtain the radio background. This is very sensitive to the evolution of 
galaxies and we shall use radio and infrared source counts to constrain the 
evolution. We start by discussing the how the observed radio flux is related 
to the luminosity in an expanding universe. 

The observed flux S u at frequency v is related to the luminosity at fre- 
quency v' — (1 + z)u by 

L v , dv' _ L„/(l + z) 
v ~ Audi dv ~ 4vr4 1 j 

where S v has units of W Hz -1 m~ 2 , L u has units of W Hz -1 , di = (l + z)do is 
the luminosity distance, and do is the physical distance at photon reception 
(m). 
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2.1 Normal galaxies 

In calculating the radio flux of normal galaxies, we shall use the 60 micron 
luminosity, L 60 , and an observed correlation between the luminosities at 1.4 
GHz, L1.4, and at 60 micron. We may then write 

r, _ {L v , / L 1A )(L lA / L m )L m (l + z) 

bv ~ 4^| • {2} 

We use the observed correlation between the luminosities at 1.4 GHz and 
at 60 micron given by Condon PS 



L 1A = 1.69 x 10^(2.58 + 1.67 a )L 60 (3) 

where both Z^o and L\ A have units W Hz -1 , and a is the spectral index at 
60 micron. Hacking et al. |29] give the distribution of spectral indices at 60 
micron which can be approximated by a gaussian distribution with mean, 
(a(L eo ,z)) given by equation 2 of P5f| , and standard deviation a = 0.5, and 
we integrate over this distribution of a in our calculation. 

To obtain the ratio (L^/L 14 ) we need to know the spectrum of nor- 
mal galaxies in the radio region from GHz down to kHz frequencies. This 
spectrum is poorly known, and so we shall model the spectrum based on 
observations of the Galaxy and other galaxies above ~ 50 MHz together 
with an estimate of the effect of free-free absorption. In the Galaxy, the 
observed spectrum is a power-law, l v oc z/~ a , with spectral index a ~ 0.9 at 
high frequencies, and a ~ 0.4 at low frequencies |3(J. If the emission from 



the Galaxy were observed from outside the Galaxy it would be modified by 
free-free absorption by the warm and hot ionized components of the inter- 
stellar medium, and by synchrotron self-absorption. Free-free absorption by 
the warm component may be expected to be patchy as based on observa- 
tions of external galaxies |3I], B2[ [£| However, this patchy absorption 
is apparently not the cause of the observed downturn of the radio spectra of 
galaxies, but rather the losses experienced by the cosmic ray electrons at low 
energies. 

The gamma-ray emission of the Galaxy demonstrates that the low energy 
spectrum of cosmic ray electrons is modified by ionization and bremsstrahlung 



losses below about 400 MeV, and cuts off below about 50 MeV 35 . Models 



can be constructed that explain both the radio as well as the gamma-emission 
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from the Galaxy. Such models then have an approximate spectrum of cosmic 
ray electrons as follows 



(E/400 MeV)- 2 - 8 > 400 MeV 

n e (E) oc { (E/400 MeV)- L8 < 400 MeV (4) 

; 50 MeV/400 MeV)" L8 < 50 MeV 

There has to be a low energy cutoff in the electron spectrum, such as exists 
in low energy protons ( ||36|| ); such a cutoff arises from the extreme losses in 
ionization and heating of the interstellar medium on the one hand, but we 
also need to note that energetic electrons can be accelerated in Supernova 
remnant shocks only at those energies for which their Larmor radius exceeds 
that of the thermal protons in a shock (f37|), and that corresponds to a 
few tens of MeV. Below this energy the electrons gain and lose energy by 



interaction with plasma waves ([p8|). 

We use the electron spectrum of Equation |], a magnetic field of 6 mi- 
crogauss and apply the standard formulae for obtaining the emission 



coefficient (W m 3 Hz 1 sr x ) for synchrotron radiation 

1 r°° 

e v = -c 3 B(wn9) n e (E)F(x)dE (5) 

2 Jo 



where x — vjv c 
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, 2 z/ g (sin#> 7 2 , (6) 

v g is the electron cyclotron frequency, (sin 9) = 0.785 for isotropic electrons, 
and C3 and the function F(x) are given in ref. f40|. The amount of absorption 
by the hot component of the interstellar medium will also depend on viewing 
angle. We assume a density of 0.01 cm -3 and a scale height of 1 kpc [f|l 



a temperature of 3 x 10 5 K and radial extent of 15 kpc giving pathlengths 
through the galaxy ranging from 2 to 30 kpc. Since the radio synchrotron 
radiation is observed to originate in approximately the same volume, we solve 
the equation of radiative transfer for the case of the emission and absorp- 
tion coefficients being independent of position throughout this volume. The 
intensity is then 

I„ = ^-[l-exp(-7i,)], (7) 

where e u is the synchrotron emission coefficient and k v and r u are the free-free 
absorption coefficient and optical depth. For a given viewing angle we model 
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the spectrum in this way neglecting synchrotron self-absorption. Assuming 
an isotropic distribution of viewing angles, we obtain the average spectrum 
in the radio region shown in Figure [TJ(a). 

Synchrotron-self-absorption will be important at very low frequencies if 
the spectrum is not cut off by other processes. Following Longair the flux 
at frequency v from a self-absorbed synchrotron source may be approximated 
by 

3 V ~ ^- 2 ^' 2 (8) 
3zV 

where m e is the electron mass and Q is the solid angle subtended by the 
source. For a source at distance d this implies a luminosity L u = A^d 2 S u at 
frequency v and solid angle Q = A pTO j/d 2 , where A pTO j is the projected area 
of the source normal to the line of sight, we obtain 

' 2/ ' iR l^W/' 2 (9) 



3^ 2 



for a self-absorbed source. For a typical galactic magnetic field of a few micro- 
gauss (6 microgauss in the solar neighborhood f39[) and typical dimensions of 
the synchrotron emitting region of normal galaxies (radius ~ 15 kpc, height 
~ 1 kpc) we find synchrotron self-absorption to be important for L\ A > 10 20 
W Hz' 1 . 

Because synchrotron self-absorption will determine the shape of the low 
frequency part of the spectrum of normal galaxies we now give a more accu- 
rate treatment. The absorption coefficient is given by 



n e (E) 
E 2 



F(x)dE, (10) 



and where the source is self-absorbed the intensity will be given by the source 
function 

e_u _ v 2 f °°n e (E)F(x)dE 

k u c 2 J °°E 2 ±[n e (E)/E 2 ]F(x)dE- 1 } 

The maximum luminosity at frequency v is then 

L u = A7rA pioj — . (12) 
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Figure 1: (a) Model radio spectrum for normal galaxies neglecting syn- 
chrotron self-absorption. Solid lines show the range of spectra for viewing 
angles 0° (top, face-on) and 90° (bottom, edge-on); dashed line shows average 
spectrum assuming isotropic distribution of viewing directions, (b) Average 
spectra for normal galaxies with luminosity L\ A = 10 19 . . . 10 25 W Hz -1 
including the effects of synchrotron self-absorption. Dashed lines show spec- 
trum of completely self-absorbed source: long dashes - accurate treatment, 
short dashes - approximate treatment. 



This is plotted as the long-dashed line in Figure |I|(b) where average spectra 
of normal galaxies are also shown for a range of luminosities. The approx- 
imate result (Equation ^) is also shown (short-dashed line) and is seen to 
give a reasonable approximation to the more accurate treatment above. In 
summary, we note that all three effects, synchrotron self-absorption, free- 
free absorption in the hot medium, and a low energy cutoff of the electron 
spectrum, contribute to cut off the spectrum at kHz to MHz frequencies. 

The 60 micron luminosity function, p(L 60 , z), is the number of sources at 
redshift z per unit co-moving volume at 60 micron luminosity Lqq per unit of 
luminosity. It has units of Mpc -3 (W Hz" 1 ) -1 . We use the local luminosity 
function, p (L m ), based on the local visibility function given by Hacking et 
al. [2~9~|. The local luminosity function is obtained from equation 3 of [2~9~| 



Po{Lm) 
Q 



2.94 x 10 28 10 Q , 

( log L m 



Y 



B z 



X 



W 



1/2 



2.5 log L 



60; 



(13) 
(14) 



where B = 1.51, W = 0.85, X = 23.96, and Y = 5.93, assuming H = 100 
km s" 1 Mpc -1 and go = 0.5 (note: in [29 Y appears to be incorrectly given 
as 6.93). 

The 60 micron luminosity function at redshift z is 



p(L 60 , z) = -j^xPo 



^60 



(15) 



f(zr\nz)j 

where g(z) and f(z) are density and luminosity evolution functions. Accord- 
ing to Gregorich et al. pUfl , g(z) and f(z) given in table 1 of Condon are 
consistent with IRAS source counts. 

The intensity (W Hz" 1 m" 2 sr" 1 ) from all sources is obtained by integra- 
tion over redshift and 60 micron luminosity: 



dV c (1 + z) f ,(L V , \ (L x 

An J dZ ^^dlJ dLe0piL ^ z) 

where dV c (m 3 ) is the element of co-moving volume. 



-60 



>(i0 



(16) 
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2.2 Radio galaxies 

The intensity due to steep-spectrum radio galaxies is obtained in a similar 
way using the radio luminosity function at 1.4 GHz 

[ ^ih^ {J ^iI dL ^ L ^(t) L -- (l7) 

We obtain the radio luminosity function from the visibility function of 
"monsters" obtained by Condon [4l| from the Auriemma et al. [44] data 



for elliptical galaxies (visibility function parameters B = 2.3, W = 0.75, 
X = 26.1, and Y = 5.47 for H = 50 km s" 1 Mpc" 1 ). Condon @ shows 



that the 1.4 GHz source counts favour strong luminosity evolution by a factor 
~ 16 at z ~ 0.8. 



For the ratio L^//L 14 we follow Falcke & Biermann [46] in assuming the 
electrons responsible for the synchrotron radiation are e ± from tc — > fi — > e 
decay, the pions being produced as secondaries in interactions of energetic 
protons with ambient matter. Energetic protons may be accelerated in radio 
galaxies to a power- law spectrum by shock acceleration |1| [18) . We make the 
usual assumption that the electrons will be in energy density equipartition 
with the magnetic field, protons, etc., and radiate in an environment of 
typically 10 to 100 microgauss (see e.g. ref. (|47|| ). The electron spectrum 
on production has a low-energy cut-off at 7 c mc 2 ~ 100 MeV (see e.g. |4S|]) 



giving rise to a break in the synchrotron spectrum at the break frequency 

3cB\ 9 

■* = < 18 > 

where 7 C ~ 200 is determined by the pion rest mass and so corresponds to the 
low-energy cut-off of the electron energy distribution. Adopting a B± = 30 
microgauss gives rise to a break at ~ 5 MHz. Below this break frequency the 
spectrum will be I/V3 (i.e. the asymptotic form of the synchrotron emission 
coefficient at low frequency for monoenergetic electrons |j4"0|). For the spec- 



trum above the break frequency we use the mean spectral index (a) = 0.75 
suggested by Condon |4"3] . This spectrum is shown in Figure 0(a). 



Because of the much higher luminosities of radio galaxies compared to 
normal galaxies, synchrotron-self absorption will be greater at low frequencies 
than in normal galaxies. Thus the contribution of radio galaxies to the 
radio background will be less at low frequencies than that of normal galaxies. 
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We may therefore use the approximate formula (Equation 0) with A pro - s ~ 
(30 kpc) 2 and B ~ 10 microgauss to obtain the maximum luminosity, and 
modify the spectrum given in Figure 0(a). The resulting spectra are given 
in Figure 0(b) for various luminosities. 



2.3 Infrared and radio source counts 

There are uncertainties in the evolution of the luminosity functions of nor- 
mal galaxies and radio galaxies, in the radio-infrared correlation used for 
normal galaxies, and in the models for the radio spectra we used. We there- 
fore decided to calculate radio and infrared source counts and compare with 
observation to ensure that the evolution of the luminosity functions used in 
calculating the radio background is consistent with the source count data. 
The number of sources with observed flux at frequency v greater than S v 

is 

dV r 



/dv r°° 
dz ~df Jl dL »'p( L »'i z } ( 19 ) 



where 



4vrd 2 



L 



S u . (20) 



+ 

Differentiating equation |19| with respect to S v we obtain 
,„ s 1 dN f , dV c d 2 T 

We use the same luminosity functions and other input we used in calcu- 
lating the radio background to obtain the radio source counts at 1.4 GHz for 
normal galaxies and radio galaxies. Resulting radio source counts of normal 
galaxies and radio galaxies are shown in Fig. 0(a) for the case of i) no evolu- 
tion and ii) evolution of the normal galaxy population according to table 1 
of Condon , and pure luminosity evolution of the radio galaxy population 
according to f(z) = (1 + z) A at all z. We compare the total source counts 
with data obtained by Condon ME . 



The calculated source counts are dominated by the radio galaxy contribu- 
tion at high S v and by normal galaxies at small S v . As can be seen, in both 
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log(i/ / GHz) 

Figure 2: (a) Model radio spectrum for radio galaxies neglecting synchrotron 
self-absorption, (b) Average spectra for radio galaxies with luminosity L\ A = 
10 22 . . . 10 29 W Hz" 1 including the effects of synchrotron self-absorption. 
Dashed line shows spectrum of completely self-absorbed source. 
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cases, the predicted source counts fall well below the data for the case of no 
evolution and are well above the data for the evolution models assumed. We 
note that the evolution of normal galaxies according to table 1 of Condon is 
also approximately pure luminosity evolution with f(z) = (1 + z) A . Our ap- 
proach to this problem is to set g(z) = 1 and modify the luminosity evolution 
to fit the source counts. We assume luminosity evolution of the form 



m = T v» ; i° (22) 




and we adjust z to fit the data. We find the best fit is given by Zq = 0.8 for 
both normal galaxies and radio galaxies. The resulting radio source counts of 
normal galaxies and radio galaxies are shown in Fig. ^|(b) and are seen to be 
in good agreement with the data. As a further check, we have calculated the 
60 micron source counts due to normal galaxies for the case of no evolution, 
evolution of the normal galaxy population according to table 1 of Condon 



, and the evolution model described above which fits best the radio source 
counts. The results are shown in Figure |] where we see that the evolution 
model we adopt gives as good a fit as table 1 of Condon B3] . 



We note that by omitting data with statistical errors larger than 20% 
in Figure HI (data with large error bars were included in the original plot 



of Gregorich et al. |pE2"f ) we see that for fluxes between 0.3 and 3 Jy both 
models predict source counts which are higher than the data. Whether this 
is an indication of a new source population contributing at S u < 0.3 Jy 
(perhaps AGN), or some systematic effect affecting the data, remains to be 
seen. However, Gregorich et al. H2J argue that their data favour the evolution 



model of Condon [^3] rather than no evolution, and our adopted evolution 
model gives as good a fit to these data as Condon's. We therefore proceed to 
calculate the mean free path for interactions of 7-rays in our calculated radio 
background for two possible cases: (a) no evolution of the normal galaxy 
population (we assume a new source population below 0.3 Jy which does not 
contribute significantly to the radio background); and (b) evolution of the 
normal galaxy population (we assume the sources with S u < 0.3 Jy are still 
normal galaxies). 
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Figure 3: Contributions of normal galaxies (dotted curves) and radio galaxies 
(dashed curves) to the total extragalactic radio source counts at 1.4 GHz 
(solid curve), (a) Lower curves are for no evolution and upper curves are for 
evolution of normal galaxies according to table 1 of Condon [4~3| and pure 
luminosity evolution of radio galaxies described by (1 + z) 4 at all z. (b) Best 
fitting pure luminosity evolution as described in the text; solid curve gives 
the total source count. Data are frorq.g'ef. 
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Figure 4: Contributions of normal galaxies to the 60 micron source counts 
for no evolution (dotted curve), evolution according to table 1 of Condon 
[fill (dashed curve), and with pure luminosity evolution as described in the 
text (solid curve). Data are points with error bars smaller than 20% from 
the summary in Figure 2 of ref. ||42|| . 
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3 The radio background and interactions of 
7- rays 

The contributions to the extragalactic radio background intensity from nor- 
mal galaxies, radio galaxies, and the cosmic microwave background, together 
with the total estimated radio background intensity are plotted in Figure [5[ 
Our result is compared with the total extragalactic radio background inten- 
sity estimated from observations by Clark et al. || and with a theoretical 
estimate made several years ago by Berezinsky ]/]]. We note that with our 
adopted model for the radio spectrum from normal galaxies, the normal 
galaxies dominate the background as suggested by their dominance in source 
counts at low flux density levels (pE3[). 

Photon-photon interactions are described in [p0| and references therein. 
The mean interaction length, A, of a photon of energy E is given by, 



[X(E)}- 1 = ^sl" de^f I ' dssa(s), (23) 

where n(e) is the differential photon number density of photons of energy 
e = hu, 

n{e) = (24) 
he hv 



and a(s) is the total cross section for photon-photon pair production |ojJ for 
a centre of momentum frame energy squared given by 

s = 2eE(l -cos 9) (25) 

where 9 is the angle between the directions of the energetic photon and the 
background photon, and 



s min = (2m e c 2 ) 2 , (26) 

e ■ - ^ C2)2 (27) 
s max (e,E) = AeE. (28) 

The interaction length for photon-photon pair production in the radio 
background is plotted in Fig. 10 along with those for competing processes 
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Figure 5: Contributions of normal galaxies (dotted curves), radio galaxies 
(long dashed curve), and the cosmic microwave background (short dashed 
curve) to the extragalactic radio background intensity (thick solid curves) 
for (a) no evolution and (b) with pure luminosity evolution as described in 
the text (upper curves), and with pure luminosity evolution only for radio 
galaxies (lower curves). Dotted band give an observational estimate of the 
total extragalactic radio background[(jntensity and the dot-dash curve 
gives an earlier theoretical estimate 0. 



and other radiation fields M. We also show the mean interaction length for 
the radio spectrum based on direct observations together with attempts at 
subtraction of the effects of galactic absorption and background 0. 

4 Conclusion 

Motivated by a new interest in electromagnetic cascades through the uni- 
verse at extremely high energies, we have made a new calculation of the ex- 
tragalactic radio background radiation down to kHz frequencies. The main 
contribution to the background is from normal galaxies and is uncertain due 
to uncertainties in their evolution. The 60 micron source counts from IRAS 
above 0.3 Jy appear consistent with no evolution provided there is a new 
source population (possibly AGN) contributing below 0.3 Jy. An alternative 
interpretation of the data is that there is strong evolution of normal galaxies 
giving agreement with the source counts above 3 Jy and below 0.1 Jy (but 
not for 0.1 < S v < 3 Jy). This gives rise to a factor of 5 uncertainty in 
the radio intensity at kHz frequencies, and this translates to a factor of 5 
uncertainty in the mean free path at 10 12 GeV. If there is a new source pop- 
ulation contributing to the infrared source counts it may also be important 
in determining the infrared background which limits the transparency of the 
universe to TeV energy gamma rays. Clearly, it is vital to determine the 
nature of the sources which dominate the 60 micron counts below 0.3 Jy. 

We calculated the radio background for the two assumptions about the 
evolution of normal galaxies, and in both cases the background we obtain 
exceeds previous estimates at low frequencies. By examining Fig. ^| we find 
that for the radio background calculated in this paper photon-photon pair 
production on the radio background is the dominant interaction process for 
photons over four or five decades of energy from 3 x 10 10 - 5 x 10 10 GeV to 
10 15 - 5 x 10 15 GeV, above which double pair production on the microwave 
background dominates. We estimate the mean free path to be ~ 1 - 5 Mpc 
at 10 12 GeV. Using the radio background estimated by Clark || photon- 
photon pair production on the radio background would only be important 
only up to 10 13 GeV, and the mean free path at 10 12 GeV would be a factor 
of 3 - 10 larger. This difference will be very important in electromagnetic 
cascades initiated by particles with energies up to the GUT scale produced 
at topological defects. 
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Figure 6: The mean interaction length for pair production for 7-rays in the 
Radio Background calculated in the present work (solid curves labelled R: 
upper curve - no evolution of normal galaxies; lower curve - pure luminosity 
evolution of normal galaxies) and in the radio background of Clark (dotted 
line). Also shown are the mean interaction length for pair production in the 
microwave background (2.7K), the infrared and optical background (IR), 
and muon pair production (/i + /i~) and double pair production (4e) in the 
microwave background [§]. 
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